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Recently there have been suggestions that the Type la supernova data can be explained using 
only general relativity and cold dark matter with no dark energy. In "Swiss cheese" models ol 
the Universe, the standard Friedmann-Robertson- Walker picture is modified by the introduction ol 
mass-compensating spherical inhomogeneities, typically described by the Lemaitre-Tolman-Bondi 
metric. If these inhomogeneities correspond to underdense cores surrounded by mass-compensating 
overdense shells, then they can modify the luminosity distance-redshift relation in a way that can 
mimic accelerated expansion. It has been argued that this effect could be large enough to explain 
the supernova data without introducing dark energy or modified gravity. We show that the large 
apparent acceleration seen in some models can be explained in terms of standard weak field gravi- 
tational lensing together with insufficient randomization of void locations. The underdense regions 
focus the light less than the homogeneous background, thus dimming supernovae in a way that can 
mimic the effects of acceleration. With insufficient randomization of the spatial location of the voids 
and of the lines of sight, coherent defocusing can lead to anomalously large demagnification effects. 
We show that a proper randomization of the voids and lines of sight reduces the effect to the point 
that it can no longer explain the supernova data. 



PACS numbers: 98.80.-k, 98.80.Es, 98.62.Py 
I. INTRODUCTION 

Supernova cosmology [l], has taught us that the Uni- 
verse is not well described by a homogeneous and flat 
matter dominated model, with gravity governed by gen- 
eral relativity. Type la supernovae at a given redshift ap- 
pear dimmer than expected based on this simple model. 
Many additions to the model have been developed to ex- 
plain this discrepancy, including dark energy and modifi- 
cations of general relativity on cosmological length scales. 
Another suggestion is that the discrepancy results from 
large scale cosmological inhomogeneities, and arises from 
an unjustified fitting of supernova data to homogeneous 
models. It is well known that inhomogeneity in the cos- 
mological matter distribution affects the redshifts and 
apparent magnitudes of supernovae [1, 0, [H, 0] , but could 
such effects be large enough to trick us into thinking that 
the Universe is accelerating when it is actually not? This 
issue has been vigorously debated in the literature; see 
Ref. for a review. 

One way to address the inhomogeneity issue in a 
fully nonlinear and relativistic fashion is with so-called 
"Swiss cheese" models [1, [tl [To|, EH- In these models, 
the homogeneous Friedmann-Robertson- Walker (FRW) 
metric is modified by the introduction of spherical re- 
gions of the spherically symmetric Lemaitre-Tolman- 
Bondi (LTB) dust spacetime. In this paper we will re- 
strict attention to such models where the FRW back- 
ground is flat and matter dominated. The spherical in- 
homogeneities must be mass compensating and comoving 
in order to obtain a self consistent solution with the LTB 
and FRW regions evolving independently. The LTB re- 



gions can be of any size and can be located anywhere, 
so long as they do not overlap. However, for a realis- 
tic model of the Universe, the sphere locations should be 
randomized, i.e. they should not placed on a lattice of 
any sort. In this paper, we show that the randomization 
issue is important and can have a dramatic effect on the 
model's predictions. 

Marra, Kolb, and Matarrese (hereafter MKM) [H, 0, 
14] use a Swiss cheese model where the LTB spheres cor- 
respond to underdense voids surrounded by overdense 
shells of matter. They place these spheres on a regular 
lattice, and compute numerically the luminosity distance- 
redshift relation along a line passing through the spheres' 
centers. They show that the relation is modified enough 
to mimic a significant dark energy density. 

We show that the results of MKM can be explained 
to ~ 10% accuracy with the formalism of standard weak 
field gravitational lensing. The special line of sight that 
they consider leads to a large dimming because this is the 
line of sight with the lowest mean density, corresponding 
to less focusing. We also show that if the void locations 
and lines of sight are properly randomized, then the mean 
effect is very small, in accordance with previous lensing 
studies [1, [TEI, [l|| [TtJ ■ Using a Monte Carlo simulation, 
we compute the statistical distribution of distance mod- 
ulus shifts at fixed redshift in the more realistic random- 
ized scenario. 

This paper is organized as follows. Section II reviews 
Swiss cheese models, focusing on the specific model con- 
structed by MKM. Section III shows that the apparent 
acceleration found by MKM is predominantly the result 
of gravitational lensing and of insufficient randomization 
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of the void locations. In Sec. IV, we show that random- 
izing the void locations reduces the effect significantly, 
to the point that this model cannot explain the apparent 
acceleration seen in the supernova data. In Sec. V, we 
provide our concluding remarks. 



II. SWISS CHEESE MODELS 

The MKM model consists of a flat, matter dominated 
FRW background with metric 



ds 2 = -dt 2 + a 2 {t) (dr 2 + r 2 dfl 2 ) 



(2.1) 



and density pFRw(t) oc 1/a 3 . From this FRW back- 
ground, spheres of constant comoving radius R are re- 
moved and replaced with spherically symmetric LTB dust 
solutions. For simplicity we describe this construction for 
a sphere centered at the origin, and for r < R. The LTB 
metric in r < R is UM 



ds 2 



-dt 2 



F /2 (r, t) 
1 + 2E(r) 



dr 2 +Y 2 (r, t)dn 2 



(2.2) 



where E(r) is the so-called energy function, primes de- 
note partial derivatives with respect to r and overdots 
denote partial derivatives with respect to t. For the case 
E(r) > considered by MKM, the angular diameter dis- 
tance Y(r,t) is given parametrically by 



Y(r,t) 



M(r) 
12irE(r) 



(cosh it — 1) 



t - t b (r) 



M(r) 

7372 



6tt [2E(r] 



(sinh u 



(2.3) 



(2.4) 



where the mass function M{r) and the bang time func- 
tion tb(r), as well as the energy function E(r), are freely 
specifiable functions. Following MKM, we use units with 
%ttG = 1. The density of the dust is given by 



M'(r) 



47rF'(r, t)Y{r,tf 



(2.5) 



There are three conditions necessary for a consistent 
matching of the LTB region with the surrounding flat 
FRW spacetime, satisfying the Israel [l!| junction condi- 
tions [2fJ. The first is that the energy function goes to 
zero at the boundary 



E(R) = 



(2.6) 



The second is that the bang time function goes to zero 
at the boundary (assuming that the big bang occurs at 
t = in the FRW background) 



t b (R) = 



(2.7) 



The third is that the LTB region be mass compensating, 
meaning that the enclosed mass has to equal the mass 



originally excised from the FRW background. This re- 
quirement is that 



M{R) 



4tv 



-PFRWfiR 



(2. 



where the constant Pfrw,o is defined by /Ofrw(^) — 
Pfrw,o/o(*) 3 - We also choose the scale factor a(t) to 
equal 1 today. 

In addition to satisfying the Israel junction conditions, 
it is convenient to specialize to models for which the coor- 
dinates (t, r) match up in a C 1 manner across the bound- 
ary r = R. This guarantees that the basis vectors d/dt 
and d/dr will be continuous across r = R, and justifies 
our use of the same notation (t, r) for the coordinates 
outside and inside the sphere. It can be shown that the 
necessary and sufficient conditions for this continuity are 

M'(R) = 3M(R)/R , t' b (R) = , E'{R) = . (2.9) 

MKM specialize to a subclass of LTB models defined 
by imposing, at some initial time t = U, the two condi- 
tions 



Y(r, U) = a(U)r , Y(r, U) = a{U)r 



(2.10) 



These conditions fix two of the three functional degrees of 
freedom of the LTB model, one of which is gauge. MKM 
parameterize the remaining, nongauge, functional degree 
of freedom in terms of the density profile at the initial 
time po(r) = p{r,ti). The mass inside radius r is given 
in terms of this density profile by 

M(r) = 4ira(U) 3 [ p (f)f 2 df . (2.11) 
Jo 

It is convenient to parameterize this as 

M(r) = M FRW (r) [1 + 5(r)} , (2.12) 

where 

4"7T 

M FRW (r) = —r 3 a{U) 3 p FRW {U) (2.13) 

is the mass enclosed that would apply in the FRW back- 
ground, and S(r) is the fractional mass excess. The en- 
ergy function E(r) and bang time function t(,(r) can be 
written in terms of the fractional mass excess as 



E(r) = 



M FRW (r) 



and 



t b (r) = U ll + 



6na(ti)r 



3[l + *(r)] 



S(r), 



(2.14) 



tanh 



28{r) 



25{r) 



(2.15) 

The last formula is valid only for S < 0, the case con- 
sidered by MKM. From the results (|2~TT|) . (jiO!?)) , and 
(I2.15[) . it follows that all of the matching and continuity 



3 



conditions (|2.6p - (|2.Q[) will be satisfied if we choose the 
fractional mass excess d(r) to satisfy 



S(R) = , S'(R) = . 



(2.16) 



The first of these conditions is the mass compensation 
condition, while the second is a requirement that the den- 
sity be continuous across r = R. 

In the next section, we will we make the approxima- 
tion of using the relativistic LTB density profile p(r, t) 
as an input to the standard weak field lensing formal- 
ism. To summarize the discussion above, the procedure 
for computing the density p(r, t) from a specified initial 
density profile po(r) is as follows, (i) Compute the mass 
enclosed M{r) using Eq. (|2.1ip and the fractional mass 
excess 5{r) using Eqs. (|2.12[) and (|2.13[) . (ii) Use the frac- 
tional mass excess to compute the energy function (j2 . 14[) 
and bang time function (|2.15p . (iii) Then compute the 
angular diameter distance Y(r, t) from Eqs. (|2.3p - (|2.4p . 
and finally insert this into the expression (12. 5p for the 
density p(r, t). 

For the initial density profile for the LTB patches, 
MKM choose the Gaussian profile 

pair) = Ap FRW (ti) {exp [-B (r - r ) 2 ] + c} (2.17) 

with constants A w 2, B « 4 x 10~ 4 Mpc~ 2 , r « 
308 Mpc, and C w 5 x 10~ 5 . This is the density at the 
initial time U — 0.2, where units are chosen such that 
the big bang occurred at t = and today is at t = 1, 
implying that the scale factor is a (t) = t 2 / 3 . Here, r 
is in Mpc and the present-day radius of an LTB patch 
is R = 350 Mpc. We will use an initial density pro- 
file that is very similar to that of MKM for all of our 
analyses below, except that we slightly modify these con- 
stants to accurately enforce mass conservation, while set- 
ting C = for simplicity. Such small changes affect our 
results only negligibly. We plot the resulting density con- 
trast A(r,t) = [p(r,t) - PFRw{t)]/pFRw{t) in Fig. [T] for 
t = 0.2, 0.6, and 1. 



III. THE RESULTS OF MARRA, KOLB, AND 
MATARRESE 



In the MKM model, the void regions are placed on 
a regular cubic lattice wherein the voids are just touch- 
ing. The chosen line of sight is throu gh the center of 
a sequence of voids; see Fig. 2 of Ref. [13| for a sketch 
of their setup. The void size is such that the distance 
through 5 voids corresponds to z « 1.8. After numeri- 
cally integrating the fully relativistic equations for light 
rays, MKM display results for the change in the redshift 
and distance modulus, both as a function of source red- 
shift, in Figs. 4 and 5 of fl3| . They find that there is 
a negligible effect on the redshift and there is a signif- 
icant effect on the distance modulus, with supernovae 
systematically becoming demagnified. This demagnifica- 
tion increases with the number of void crossings. The 
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FIG. 1: The density contrast A(r,t) = [p(r,t) - 
PFRw(£)]/pFRw(i) as a function of LTB coordinate r in the 
model of MKM at the times t = 0.2 (blue, bottom curve), 0.6 
(pink, middle curve), and 1 (yellow, top curve). 



effect is large enough that supernovae that are 5 voids 
away suffer a distance modulus shift of about 0.34 mag, 
and this suffices to move them to the Hubble diagram of 
a J7m = 0.6, f^A = 0.4 universe. 

We now argue that this result is predominantly due 
to gravitational lensing. Although the inhomogeneities 
are mass compensating, light beams passing through 
the spheres' centers encounter less mass than beams 
that intersect the spheres with larger impact parame- 
ters. Therefore, centered beams are focused less than 
they would be in the FRW background. 

In weak lensing theory, the net magnification p pro- 
duced by a density perturbation Ap(x, 77) is p, = (1 — 
k)~ 2 . Note that we are neglecting shear, which gives a 
contribution to p of the same order as the k 2 corrections. 
The lensing con verg ence k is given by the integral along 
the line of sight [1l| 



V ° j (V-Vs)(VO ~V) f n2a r / n 1 

a(V) Ap[x(?7),77] . 



k = AttG / dri 

>vs 



VO - TjS 



Here the integral is with respect to conformal time 77 and 
is along the unperturbed ray x(ry), and rjs and rjo are 
the values of conformal time at the source and at the ob- 
server. Rewriting this integral in terms of the comoving 
distance x = rjo — tj, expressing it in terms of the den- 
sity contrast A(x, 77) = Ap/ppRw, and eliminating pfrw 
using the Friedmann equation gives 



„ , ' S dx X{ * S X \ A(x(x), r)o-x), (3.2) 




where x$ = fjs ~ Vo- Next, since the individual spheres 
are small compared with the horizon scale, we can treat 
the first factor in the integrand as constant across each 
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sphere, and we can also neglect the time dependence of 
A within each sphere to leading order. This gives 



K= 3 2 y x^xs-x ) 
2 x s a( Vj ) J 



with 



dx A(x(x),^) , 



(3.3) 



(3.4) 



where Xj is the value of x at the point of closest approach 
of the ray to the center of sphere j, and rjj = rjo — Xj 
is the conformal time at that point of closest approach. 
The interval Sj is the intersection of the ray with the 

jth sphere, given by \x — Xj\ < \J R 2 — b 2 , where bj is an 

impact parameter. Rewriting the integral (|3.4[) in terms 

of the comoving radial coordinate f = ^Jb 2 + (x — Xj) 2 

gives 



Xj = 2 / df 



2 -b 2 



A(r,7/j-) 



(3.5) 



To evaluate this integral, we use the relation f = 
Y(r,tj)/a(tj) between the comoving radial coordinate f 
and the LTB radial coordinate r, where tj = t(r)j). Fi- 
nally the distance modulus shift Am is related to the 
lensing convergence k by 



Am = 51og 10 (1 — ft) 



(3.6) 



We have computed the distance modulus shifts Am 
using Eqs. (O and for the model of MKM, with 

our aforementioned slight modifications. We also use 
Ha = 70 km/s/Mpc. These results are summarized in 
Table HI Comparing with Fig. 5 of Ref. [HI, we see that 
our weak field results reproduce their fully relativistic re- 
sults to ~ 10% accuracy. This is consistent with what 
one would expect given the approximations that we have 
made, since RH„ ~ 0.1 for R = 350 Mpc and H = 70 
km/s/Mpc. We conclude that gravitational lensing pro- 
duced by the coherent arrangement of voids provides a 
simple physical explanation for the phenomenon seen in 
their analysis. 



TABLE I: Supernova distance, in terms of the number of voids 
in between source and observer, and our result for the distance 
modulus shift Am if the voids are in a lattice. The void radius 
is R = 350 Mpc. 



Distance (in voids) 


Am 


1 


0.015 


2 


0.047 


3 


0.11 


4 


0.20 


5 


0.33 



We can get a good estimate the coherent MKM effect 
by using the Dyer and Roeder [22], [23[ formalism. This 



is based the work of Zel'dovich [24|, wherein he consid- 
ered cosmological observations in a flat matter dominated 
universe where there was no matter inside any observed 
light beams, and hence no gravitational focusing. Ignor- 
ing the effects of shear, the observed flux from a source 



at redshift z in such a model is F, 



empty 5 



where 



F 



empty 



25[l-(l + z)- 1 < / 2]2 



FfUled (l + z) 2 [l-(l+z)-5/ 2 ] 2 



(3.7) 



with F fui e d being the flux for a filled beam. Thus, high 
redshift sources would be dimmed. The corresponding 
magnitude shift is 



Am 



empty 



5 log 



10 



(l + z)[l-(l + z)- 5 / 2 ] 
5[1 - (1 + z)- 1 / 2 ] 



(3.8) 



and for a model in which we observe through a complex 
of voids, this must be an upper bound to the magnitude 



shifts possible. For z = 1.8, we find Am, 



empty 



0.546. In 



the absence of shear, no arrangement of matter along the 
line of sight can produce a magnitude shift larger than 
this. 

Dyer and Roeder generalized this analysis to partially 
filled beams. They found 



Am Q = 5 log 



in 



(1 + z^-^l^l - (1 + z)-?/ 2 

P[l - (1 + Z)-V2] 



(3-9) 

where j3 = V25 — 24a and a is the filling factor of the 
beam: a — 1 corresponds to a flat dark matter dominated 
model and a = corresponds to empty beams. 

The limiting case is a — for empty beams, but what is 
the limiting case for passage through a mass compensated 
void? Consider an extreme example in which the void 
has all of its mass concentrated in a shell at its boundary. 
Then if the mass is M and the radius is i?, the integrated 
column density through the sphere is M/2irR 2 = 
i.e. this is zero except for the shell, which is pierced twice. 
This is 1/3 of the value for a smooth sphere, and hence 
for this model a — 1/3 would be appropriate. Using this 
value we get (3 — vT7 and 



ATOi/3 = 5 log 



10 



(l + z )(Vl7-l)/4[ 1 _( 1 + z )-^7/2] 



17[1- (1 + z)- 1 / 2 } J 
(3.10) 

For z = 1.8, this gives Am^ = 0.369. This just slightly 
overestimates the MKM result of Am w 0.34. 



IV. VOID RANDOMIZATION 

As we have argued, the coherent arrangement of the 
voids on a lattice in the model of MKM artificially boosts 
the demagnification effect. In this section, we show that 
the net effect of the voids is small if the void locations 
relative to the line of sight are sufficiently randomized. 

As is well known, the mean magnification due to gravi- 
tational lensing should be negligible on average in models 
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such as ours with small, mass-compensatin g in homogen- 
ities that are randomly distributed d, [TEI. Il7j. We can 
understand this result in a different way with a simple ar- 
gument. First, consider a spherically-symmetric matter 
distribution with a total radius R, a total mass M, and 
a radius-dependent density p(r). Let us further assume 
that a small beam comes in from the z direction and then 
encounters the sphere with an impact parameter b, which 
will then be in the x — y plane. The magnification of the 
beam is proportional to the integrated column density 
along the beam's unperturbed path, which is defined to 
be 



E = 



p[r{\)] dX 



(4.1) 



path 



where the path depends on 6, and A is the affine param- 
eter of a central ray. Since we are using the unperturbed 
path for this, it is clear that we do the above integral 
along a straight line, parallel to the z axis. Then the 
total mass of the spherical distribution comes from inte- 
grating E over the remaining two spatial axes, the x and 
the y: 

M = / p(r)dV = / p(r)dxdydz 

J all space J all space 



E [b(x, y)] dxdy 



x — y plane 
R 



f 

/ E(r)27rr<ir . 
Jo 



(4.2) 



The average focusing resulting from a spherical mass dis- 
tribution is proportional to the average over impact pa- 
rameters of the integrated column density, 



(E) = 



1 

R 2 



E(6) x 2bdb 



1 

R 2 



T,(r)2rdr 



1 (M 



M 



(4.3) 



which only depends on total mass and size of the spheri- 
cal region, and not on how the mass is distributed. Thus, 
we see that removing a sphere of FRW spacetime and re- 
placing it with a mass-compensating LTB patch, as we 
did in the model discussed above, will not affect the mag- 
nification on average. This is true within the domain of 
validity of weak lensing theory. The distance modulus 
shift depends on the logarithm of the magnification and 
thus it is not necessarily going to have a vanishing mean, 
but we nonetheless find that it is small in the model con- 
sidered here. 

To find the statistical distribution of the distance mod- 
ulus shifts that one would expect after passage through 5 
voids in a randomized version of this Swiss cheese model, 
we do a Monte Carlo analysis. In this analysis, we pick 
random impact parameters for each of the 5 void pas- 
sages, with b 2 being uniformly distributed between and 
R 2 . We find that the statistical distribution of distance 




FIG. 2: A histogram of the distance modulus shifts Am found 
for 10, 000 realizations of our randomized void scenario 



modulus shifts has a mean of —0.003 mag and a stan- 
dard deviation of 0.11 mag. Note that we have neglected 
shear, which might change the mean magnitude shift by 
a factor of order 2. Figure [5] shows the distribution of 
these shifts found for 10, 000 realizations of voids with 
random impact parameters relative to the line of sight to 
the light source. 



V. DISCUSSION 

MKM found significant changes to the distance mod- 
ulus of distant supernovae in a particular Swiss cheese 
model. We have argued that the large size of their ef- 
fect stems from the fact that their voids are on a regular 
lattice, with the line connecting observer and source pass- 
ing straight through the void centers. Their result can be 
reproduced with ~ 10% accuracy using weak field gravi- 
tational lensing, and can be explained by noting that the 
integrated column density is lowest along lines of sight 
that pass through the centers of the voids. The lower 
column density leads to less focusing and thus dimmer 
supernovae. 

We also showed that in a more realistic Swiss cheese 
model, the mean magnification due to gravitational lens- 
ing is very small. We simulated the distribution of dis- 
tance modulus shifts obtained from a model where the 
light rays no longer pass only through the centers of the 
voids, but instead encounter them with random impact 
parameters. For N — 5 voids out to z — 1.8, MKM 
obtained a coherent distance modulus shift of 0.34 mag, 
whereas our randomized model gives a mean of —0.003 
mag and a standard deviation of 0.11 mag. This mean 
is small and this standard deviation is comparable to the 
intrinsic scatter of Type la supernovae, and thus this 
particular model cannot explain the supernova data. An 
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analysis of the luminosity distance-redshift relation in a 
scenario with randomly-sized voids is the subject of fu- 
ture work. 

While our results clarify the significance of the special 
setup underlying the specific model of MKM, they do not 
exclude the possibility that inhomogeneties may account 
for some of the apparent acceleration of the Universe. 
However, our work does suggest that coherent under- 
dense structures on scales of order the horizon size gen- 
erally are needed to obtain large apparent accelerations. 
Examples of models with such very large-scale structures 
arc the spherical void models of Refs. ^1 [H, [H El, EH . 
A general argument against such models is fine tuning: 
one would generically expect large anisotropies in the lu- 
minosity distance at a given redshift in any model with 
a large apparent acceleration, and the observational lim- 
its on anisotropy are small. Going beyond fine tuning, 
there are limits on large spherical void models from cur- 
rent observations of the cosmic microwave background 
spectrum [3(| and low redshift Type la supernovae [3l| . 
and further constraints could come from future observa- 
tions of baryon acoustic oscillations [32] and time drifts 
of cosmological redshifts [HI, [34[ ■ 

We also note that in spherical void models, fitting lumi- 



nosity distance data to spatially flat cosmological models 
would generally suggest the presence of dark energy. If 
we allow more general fits to non flat models, then we 
will get a fit with non zero curvature and non zero dark 
energy, with the relative proportion depending on the de- 
tails of the model. In linear theory, the shrinking mode 
of an LTB model is the same as the mode parameter- 
ized by the bang time function, and the growing mode 
is the mode parameterized by the energy function, which 
parameterizes curvature. Therefore in a model with no 
important shrinking mode component, it is likely that 
the best fit will be have a higher proportion of curvature 
than dark energy. 
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